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■ Abstract 

We calculate the nucleon parameters in isospin symmetric and asym- 
metric nuclear matter using the QCD sum rules. The higher moments 
of the nucleon structure functions are included. The complete set of the 
nucleon expectation values of the four-quark operators is employed. We 
analyze the role of the lowest order radiative corrections beyond the loga- 
rithmic approximation. 
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1 Introduction 

We investigate the vector and scalar self-energies of nucleons in nuclear matter 
composed by the neutrons and protons, distributed with densities p n and p p . We 
calculate the dependence on the total density p = p p + p n and on the asymmetry 
parameter /3 = (p n - p p )/ (p p + p n ). 

The QCD sum rules were invented in paper pQ to express the hadron param- 
eters through the vacuum expectation values of QCD operators. Being initially 
used for the mesons, the method was expanded in [2] to the description of the 
baryons. The approach succeeded in describing the static characteristics as well 
as some of the dynamical characteristics of the hadrons in vacuum — see, e.g. 
the reviews [3]. 

Later the QCD sum rules were applied for investigation of modified nucleon 
parameters in the symmetric nuclear matter [1]. They were based on the Borel- 
transformed dispersion relation for the function il m (g) describing the propagation 
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of the system with the quantum numbers of the nucleon (the proton) in the 
nuclear matter. Considering nuclear matter as a system of A nucleons with 
momenta pi, one introduces the vector p = (EpA/A, which is thus p ~ (m, 0) in 
the rest frame of the matter. The function H m (q) can be presented as H m (q) = 
n m (g 2 , (pip, q)) with the arbitrary function (p(p,q) being kept constant in the 
dispersion relations in q 2 . 

The general form of the function U m can thus be presented as 

n m (g) = q^Jq 2 , s) + p^W m {q\ s) + IU^q 2 , a). (1) 

The in-medium QCD sum rules are the Borel-transformed dispersion relations for 
the components lP m iq 2 , s) (j = q,p, I) 

The spectrum of the function IT m (g) is much more complicated than that of 
the function il (g) describing the propagation of the system with the quantum 
numbers of the nucleon in the vacuum. The choice of the function (pip, q) is 
dictated by the attempts to separate the singularities connected with the nucleon 
in the matter from those connected with the properties of the matter itself. Since 
the latter manifest themselves as singularities in the variable s = ip + q) 2 , the 
separation can be done by putting (pip, q) = ip + q) 2 and by fixing [5] pip, q) = 
ip + q) 2 = s — 4m 2 im is the nucleon mass). 

By using Eq. ([2]) the characteristics of the nucleon in nuclear matter can be 
expressed through the in-medium values of QCD condensate. The possibility of 
extension of "pole + continuum" model [H [2] to the case of finite densities was 
shown in [5] [7]. 

The lowest order of OPE of the lhs of Eq. (TS]) can be presented in terms of the 
vector and scalar condensates [5], [7J. Vector condensates v % = {M\q l ^^q l \M) 
of the quarks with the flavor i (|M) denotes the ground state of the matter) are 
the linear functions of the nucleon densities p n and p p . In the asymmetric matter 
both SU(2) symmetric and asymmetric condensates 

v, = (M |u(0) 7m u(0) + d(0)^d(0)\M) =v; + v1 

«(-> = (M|tt(0) 7 „ti(°) - <*(0) 7 „d(0) |M> = «« - «J* 
obtain nonzero values. In the rest frame of the matter v l = v 1 5^q, = vS^o, 
v u~ = v ^^o- We can present v % = (n\q l 'joq t \n)p n + {p\q l ^oq l \p)p p - The val- 
ues (iV|(f 7 g*|iV) are just the numbers of the valence quarks in the nucleons 
(nlujouln) = (p|g?7og?|p) = 1, (p|«7ow|p) = {n\d^od\n) = 2, and thus 

v u = p n + 2p p = p ( I - | j , v d = 2 Pn + Pp = p(l + %). (3) 
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Hence, we obtain v(p) = vnP, \p, (3) = f3v N p with — 3, v^' — — 1. 

The lhs of Eq. (JSJ) contains the SU(2) symmetric scalar condensate n m (p) = 
(M\u(0)u{0)+d{0)d(0)\M) , and SU(2) asymmetric one ( m {p,P) = (M\u(0)u(0)- 
d(0)d(0)\M) . These condensates can be presented as 

K>m(p) — «o + k(p) j n(p)=K N p + ... , k n = (N\uu + dd\N) , (4) 

= ^m(O) is the vacuum value, and 

Cm(p, P) = ~P(CnP + ■■■), Cn = (p\uu - dd\p) . (5) 

The dots in the rhs of Eqs. (jl]) and ([3]) denote the terms, which are nonlinear in 
p. In the gas approximation such terms should be omitted. The SU(2) invariance 
of vacuum was assumed in Eq. ([5]). The expectation value is related to the 
ttN sigma term u^n [H]- 
The gluon condensate 

g m (p) = (M\—G 2 (0)\M) = g + g(p) = g + g N p+ ... 

7T 

9o — fl'm(O) is the vacuum value and g^ = — | m obtained in a model-independent 
way. 

We shall analyze the sum rules in the gas approximation. It is a reasonable 
starting point, since the nonlinear contributions to the most important scalar 
condensate n(p) are relatively small at the densities of the order of the phe- 
nomenological saturation density po — 0.17 fm -3 of the symmetric matter [7]. 

In the second part of the talk we discuss the role of the radiative corrections. 
The analysis [2] included also the most important radiative corrections, in which 
the coupling constant a s is enhanced by the "large logarithm" lng 2 . The correc- 
tions (a s In q 2 ) n have been included in to all orders for the leading OPE terms. 
This approach provided us with good results for the nucleon mass and for the 
other characteristics of nucleons. 

However, inclusion of the lowest order radiative corrections beyond the log- 
arithmic approximation made the situation somewhat more complicated. A nu- 
merically large coefficient of the lowest radiative correction to the leading OPE 
of the polarization operator n (g) was obtained in [9]. A more consistent calcu- 
lation [PJ provided this coefficient to be about 6. Thus, the radiative correction 
increases this term by about 50% at \q 2 \ ~ 1 GeV 2 , which are actual for the SR 
analysis. This uncomfortably large correction is often claimed as the most weak 
point of the SR approach [9j. 

The radiative corrections of the order a s In q 2 and a s for the contributions up 
to q~ 2 have been calculated in [TO] . 

The further development of the nucleon SR in nuclear matter needs the cal- 
culation of the radiative corrections. This work is in progress and now we have 
present the analysis of the role of the radiative corrections in vacuum [11] . 
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2 Sum rules in nuclear matter 



We present the nucleoli propagator in nuclear matter as 

G], 1 = q^-m-Z, E = q^Zg+p^^ + Zs (6) 

with the total self-energy E . We shall use the QCD sum rules for the calculation 
of the nucleon characteristics 

v> * 777, + E s 

E„ = — — , m = - _ - , E s = m - m , (7) 

identified with the vector self-energy, Dirac effective mass and the effective scalar 
self-energy — see e.g. [T2] . 
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q 2 -w? m 

with E„ and m* defined by Eq. (JTj). The new position of the nucleon pole is 

m 1 + K/m ' N (l-E 9 )(l + E„/m) ' { > 

We present also the result for the single-particle potential energies 

U = E: + £„ . (10) 

We trace the dependence of these characteristics on the total density p and on 
the asymmetry parameter. 

The Borel-transformed sum rules take the form 

L^M^Wl) = A m (M 2 ), (11) 
W m {M\W 2 m ) = -S,A m (M 2 ), (12) 
LUM^Wl) = m*A m (M 2 ) (13) 

with A m (M 2 ) = A^e - ™" 1 ^ 2 , where A^ 2 is the effective value of the nucleon 
residue in nuclear matter, M is Borel mass (0.8 < M 2 < 1.4GeV 2 ). 

We shall include subsequently the contributions of three types. The terms 
£ J m (M 2 ) (J = p,q,I) stand for the lowest order local condensates. These con- 
tributions are similar to simple exchanges by isovector vector and scalar mesons 
between the nucleon and the nucleons of the matter. The terms u J m (M 2 ) are 
caused by the nonlocalities of the vector condensate. They correspond to the 
account of the form factors in the vertices between the isovector mesons couple 
to the nucleons. Finally, ^(M 2 ) describes the contributions of the four-quark 
condensates. They correspond to the two-meson exchanges (or to exchanges by 
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four-quark mesons, if there are any) and to somewhat more complicated structure 
of the meson-nucleon vertices. Thus we present the lhs of Eqs. (jlip - (IT3"|) as 

LL = 4 + < + <• (14) 
Actually, we shall solve the sum rules equations, subtracting the vacuum effects 

V = f + u 3 + w 3 (15) 

with £ 3 = £ 3 m — f , u 3 = u 3 m , us 3 = ijj 3 m — lUq, while £ 3 and Uq are the corresponding 
contributions in the vacuum case. 



A. Local condensates of the lowest dimensions. 

The terms t have the form: 

P = f«(M 2 ,Wl)v«(p) + f« g (M 2 ,Wl)g(p), 
P = g(M 2 } W^)v"(p,(3) , 
f = fl{M\Wl)t\p,{3), 

with the dependence on p and j3 being contained in the factors 



(16) 
(17) 
(18) 



v«(p) = 3p, v*{p,l3)=?>p(l-^\, t I (p,(3) = p(K N + C N (3) } (19) 

and the function g(p), given by Eqs. (jlj) and (J5]). The other functions are [1 

8tt 2 (s - m 2 )M 2 E 0m - M A E lm 



tt{M\W 2 m , 



m 



L 4 /9 



n 2 M 2 E 



Qui 



£4/9 ' 

8tt 2 4M 4 E ir 



3 L 4 / 9 
-4ir 2 M A E lm . 



(20) 



k=0 n 



The notation E^ m (k = 0, 1) means that the functions E(x) = 1 — e X H 
depend on the ratio W^/M 2 . The factor L = ln(M 2 /A 2 )/ ln(z/ 2 /A 2 ), where 
A = 0.15 GeV, v = 0.5 GeV. 

B. Inclusion of the nonlocal condensates. 

Finally, the higher moments and higher twists of the nucleon structure functions 
provide the contributions u % to the L l of the sum rules — Eq. ( I15p 



u q {M 2 



(u%AM 2 )+(3u%JM 2 ))p 



u%JM 2 



57T 



3L 4 / 9 



rn 



~m 2 M 2 E 0m (ria) + ^m 2 {s - m 2 )(£) 
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u%JM 2 



8tt 2 



3L 4 /9 



m 



m 2 M 2 E Gm {{rfa)-{r, d a)) 



(21) 



u p {M 2 ) = 

u p NA (M 2 
12 



K 5l (M 2 ) + (3u p N ^M 2 ))p 

r 2 



57T 



3L 4 /9 



5 (M 4 £ lm - (s - m 2 )M 2 E 0m ) {w) 
18 



^-m 2 M 2 E 0m (va 2 ) + —m 2 M 2 E 0m (O 



u p N ^{M 2 
9 



57T 



3L 4 / 9 



3 (M 4 E lm - (s - m 2 )M 2 E 0m ) (( V u - rj d )a) 
27 



+ ^m 2 M 2 E 0m ((r ] u - V d )a 2 ) - ^-m 2 M 2 E 0m ((C ~ O) 
5 10 



(221 



(23) 



Here we denote = -0.24, £ d = -0.09 P5] and £ = £ u +£ d . We use the structure 
functions rj, obtained in [TH] for the calculation of the terms u q and u p . Here we 
denote (/) = J daf(a). 



C. Inclusion of the four-quark condensates. 

The calculations of the contributions of the four-quark condensates require the 
model assumption on the structure of the nucleon. The complete set of the four- 
quark condensate was obtained in [T7j by using the perturbative chiral quark 
model (PCQM). There are three types of contributions to the four-quark con- 
densate in the framework of this approach. All four operators can act on the 
valence quarks. Also, four operators can act on the pion. There is a possibility 
that two of the operators act on the valence quarks while the other two act on 
the pions. Using the complete set of the nucleon four-quark expectation values 
[T7| . we obtain 



(n) 49 = Ul q {P)\ + A%{p)^l 



with a = — (27t) 2 (0|mu|0). The calculations give 



^Aq 



0.11-0.21/3, A\ q = -0.57 + 0.09/3, A{ a = 1.90 - 0.92/? 



(24) 



(25) 



The contributions of the four-quark condensates to the lhs of the Borel trans- 
formed sum rules (fT5|) can be presented as 



u 3 N P 



< = 4Mfi 



Uq ' 



4 g — — 07T a, 



m 



iq 



j7T 



2m 



fi 



Aq 



57T ma. 



(26) 
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Figure 1: Results for m* and U(p) both in symmetric and neutron matter. Solid 
curves are for symmetric matter; proton (dashed) and neutron (dotted) curves 
are for the neutron matter. 



In Fig. 1 we show results for effective nucleon mass m* (J7|) and for the single- 
particle potential energies ( ITUl) in symmetric and asymmetric nuclear matter. 
During calculation we use the value kn = 8 (jl]), and (n = 0.54 (jSJ) in Eq.(19) 

PI- 



3 The role of radiative corrections 

We have made [TT] the analysis of radiative corrections to the nucleon SR in 
vacuum. The lowest OPE terms of the operators n^g 2 ) and U^q 2 ) in vacuum 
can be presented as 



Y\i 



A + A A + A e + A 8 



n 



B 3 + B 7 + Bq 



(27) 



Here the lower indices show the dimensions of the condensates, contained in the 
corresponding terms, A is the contribution of the free quark loop. 

We consider the radiative corrections to the terms A , Aq, and B 3 . The 
radiative corrections to the other terms are not included since the values of the 
corresponding condensates are known with poor accuracy. We find, following [18J, 
for the corresponding contributions from a s -corrections to the Borel transformed 
SR: 



A (M 2 ,W 2 ) = M 6 E 2 



7T V 12 



TV 



M 4 W 2 [ 1 + 



3W r 
AM 2 



-W 2 /M 2 



+ M 6 £ (-W 2 /M^ 
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A 6 (M 2 ,W 2 



3° 



(28) 



x 



l-—\- + -\hi- 



W 2 



71 \6 
B 3 (M 2 } W 2 ) 



V 



£(-W 2 /M 2 



2(/.U'£', ( 1 + - — 

2 7T 



with £(x) = J2n=i xTl /{ n ' n 0- Some terms in (28) differ from those in [18]. The 
numerical difference is, however, not very important. 

The parts of equations proportional to a s are the a s -corrections to main terms. 
Now we compare the nucleon parameters obtained as solutions of nucleon SR 
without a s -corrections, 

m = 0.930 GeV, A 2 = 1.79 GeV 6 , W 2 = 2.00 GeV 2 , 

and with a s -corrections, 

m = 0.94 GeV , A 2 = 2.00 GeV 6 , W 2 = 1.90 GeV 2 . 

These results are presented for a s =0.35. We show that in vacuum the radiative 
corrections modify mainly the values of the nucleon residue, while that of the 
nucleon mass suffers minor changes. 
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